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. The 4-th order Runge-Kutta method in the complex plane is proposed 

for numerically advancing the solutions of a system of first order differential 
equations in one external invariant satisfied by the master integrals related 
CO ■ to a Feynman graph. Some results obtained for the 2-loop self-mass MI are 

I reviewed. The method offers a reliable and robust approach to the direct 

and precise numerical evaluation of master integrals. 
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1. Introduction 

The aim to more precise and handy calculations of radiative corrections 
push for a restless research and development of new methods. 

Today the organization of the calculations is usually based on the inte- 
gration by part identities and on the evaluation of the master integrals (MI) 

In this frame the differential equations for the MI, or Master Differential 
Equations (MDE), can be used not only for their analytic calculations when 
the number of the variables and parameters is small, but also for their di- 
rect numerical evaluation, when the large number prevents the success of an 
analytic calculation, in alternative to the more commonly used integration 
methods or to the more recently introduced difference equations method. 
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Here is presented a method to get a numerical solution of the MI for 
any values of the variables and parameters, using directly the MDE and the 
values of the MI known in a given set of the variables and parameters. Some 
results obtained for the 2-loop self- mass MI are reviewed. 



Starting from the integral representation of the Nmi MI, related to a 
certain Feynman graph, by derivation with respect to one of the internal 
masses rrii [2] or one of the external invariants Se and with the repeated 
use of the integration by part identities, a system of Nmi independent first 
order partial MDE is obtained for the Nmi MI. For any of the Sg, say Sj, 
the equations have in general the form 



where Fk{n,mf,Se) are the MI, Kk{mf,Se) and Mk^i{n,mf, Se) are polyno- 
mials, while Tk{n, m?, Se) are polynomials times simpler MI of the subgraphs 
of the considered graph. The roots of the equations 



identify the special points, where numerical calculations are troublesome. 
Fortunately analytic calculations at those points come out to be possible in 
all the attempted cases so far. They might not be simple and often require 
some external knowledge, like the assumption of regularity of the solution 
at that special point. 

To solve the system of equations it is necessary to know the MI for a 
chosen value of the differential variable, sj in Eq.(^. For that purpose 
we use the analytic expressions at the special points, taken as the starting 
points of the advancing solution path. Moreover starting from one special 
point, not only the values of the MI are necessary, but also their first order 
derivatives at that point. That is because some of the coefficients Kk(mj, Sg) 
of the MI derivatives in the differential equations Eq.l^Q) vanish at that 
point. Therefore also the analytic expressions for the first derivatives of MI 
at special points are obtained, but this usually comes out to be a simpler 
task (unless poles in the limit of the number of dimensions n going to 4 are 
present). 

Enlarging the number of loops and legs increases the number of param- 
eters, MI and equations, but does not change or spoil the method. 



2. Master Differential Equations 




k,l = 1, ...,Nmi 



(1) 




(2) 



ustron printed on February 1, 2008 



3 



3. The 4-th order Runge-Kutta method 

Many methods are available for obtaining the numerical solutions of a 
first-order differential equation (4| 

The 4-th order Runge-Kutta method is a rather precise and robust approach 
to advance the solution from a point x„, where the solution y„ is known, 
to the point Xn+i = 2;„ -|- A By suitably choosing the intermediate points 
where calculating f(x,y) one obtains the 4-th order Runge-Kutta formula 

ki = Af{xn,yn), 

A ki 

k2 = Af{Xn + —,yn + —), 

k3 = A/(x„ y,y„ + y), 

ki = Af{Xn + A,yn + k3), 

ki k2 ks A;4 , . t:, , 

D 3 3 D 

which omits terms of order A^. 

To avoid numerical problems due to the presence of special points on the 
real axis, it is convenient to choose a path for advancing the solution in the 
complex plane of x. 

The extension from one first-order differential equation to a system of 
Nmi first-order MDE for the Nmi MI is straightforward 1^. 

4. Results 

I report here some results obtained with the outlined method. 

The simplest nontrivial application of the method is to the 2-loop self- 
mass with arbitrary internal masses, which has three 2-loop topologies 
the sunrise, shown in Fig^ with 4 MI, the one with 4 propagators, shown 
in Figl^l with one more MI, and the one with 5 propagators, shown in 
Fig|21 with one even more MI. The simpler 2-loop self-mass amplitudes of 
Fig |ll2l^^ can be written in the integral form as 

A{n, 'mf,p^, -ai, -02, -03, -"4, -as) = A**^"^" 
d^'ki f <Fk2 1 



D, = {kf + ml), i = 1,2,3,4,5. (5) 
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Fig. 1. The general massive 2-loop sunrise self-mass diagram. 
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Fig. 2. The general massive 2-loop 4-denominator self-mass diagram. 

where ki is the 4-momentum of the line of mass mj and 4-momentum con- 
servation is understood in the vertexes. The arbitrary mass scale /x accounts 
for the continuous value of the dimensions n; as one of the natural scales 
of the problem is the 3-body threshold of the sunrise amplitudes, we take 
usually = mi + m2 + m^. 

The 4 MI of the sunrise, the one of the 4-denominator and the one of 
the 5-denominator can be chosen as 
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the other amplitudes are related to these by integration by parts identities 



Fig. 3. The general massive 2-loop 5-denominator self-mass diagram. 
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OinilZI, which are used also to obtain the system of MDE for the MI in the 
form of Eq. (QJ . 

In the MDE of the sunrise MI the only lower order MI entering is the 
one related to the 1-loop vacuum graph 

T{n, m )- J ^2^^„_2 + ^2 " _ 2)(n - 4) ' ^' 

while for the 4-denominator MDE is necessary also the knowledge of the 
1-loop self-mass and the 2-loop vacuum MI, known analytically, and the 
sunrise MI. The 5-denominator MDE requires also the knowledge of the 
4-denominator MI. 

The function C(n) = {2^/^^)^'^ E (3 — which appears in the expres- 
sions for the MI as an overall factor with an exponent equal to the number 
of loops, is usually kept unexpanded in the limit n ^ 4, and only at the 
very end of the calculation for finite quantities is set C(4) = 1. 

When the MI are expanded in (n — 4), for k = 0,1,2,3,4,5 and i = 
1,2,3,4,5 



+ ^^/t'\^lp') + 4'\rnlp') + 0{n - 4)} , (8) 



the coefficients of the poles can be easily obtained analytically for arbitrary 
values of the external squared momentum p^, 

Fr'imlp-) = -limj+ml + ml), 

Ff^>{mt,p') = \\^^ + \{ml + ml + ml) 

-fl + mi log! -I 
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1 1 , I m 
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Fk '{rnty) = s , fc = l,2,3 



Fr^'imi.p') = +1 



Ft\mlp') = -1 - i5W(m?,ml,/) , (9) 
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where S^^\mi, ml, p'^) is the finite part of the (n— 4) expansion of the 1-loop 
self-mass 

5(n,m?,mi,p2)=C(n)|-i^-l^+5W(m?,m2,p2) + o(„_4)| , (iq) 

which is known analytically. The 5-denominator is not divergent for n ~ 4, 
so F^~'^\mf,p^) = F^~^\mf,p'^) = and F^{n,mf,p^) = F^^\inf ,p'^). 

The finite parts Flf\mf,p^) satisfy MDE of the type of Eq.©. 

The special points are easily obtained from the coefficient Kk{mf ,p'^). 
For the sunrise MI they are at = 0, oo and the roots of 



D{mj,p'^) = + (mi + m2 + ms)^ p^ + {mi + m2 - msY 
p^ + (mi — m2 + m-s)^ p^ + (m-i — m2 — m^Y 



, (11) 



For the sunrise MI the analytic expressions for their first order expansion 
were completed around the special points j3 |H1 El = 0; = oo; 

p^ = —(mi + m2 + ms)^, the threshold; p'^ = —(mi + m2 — m^)"^, the 
pseudo-thresholds. For the 4-denominator are special points also the 2-body 
threshold p'^ = —(mi -|- m^y and pseudo-threshold = —(mi — m^)'^; the 
expansion at = is completed in and \12i ll'6\ . The 5-denominator 
has special points at p'^ = 0, p'^ = oo and in some mass combinations of 
2-body and 3-body thresholds and pseudo-thresholds. 

To obtain numerical results for arbitrary values oip^, a 4th-order Runge- 
Kutta formula is implemented in a computer code, with a solution advancing 
path starting from the special points, so that also the first term in the 
expansion is necessary. 

The path followed starts from = and moves in the lower half complex 
plane of = p^/^j?, to avoid proximity to the other special points, which 
can cause loss in precision. Values between special points can be safely 
reached through a complex path. Results for arbitrary values of p^ and of 
the masses are obtained with this method for the sunrise MI in [71 llUj and 
jllj , for the 4-denominator MI in and |121 113j , for the 5-denominator 
MI in HU. 

For values oi pi very close to a special point (~ 10~^) the method fails. 
In the case of the sunrise MI we start the Runge-Kutta path from the ana- 
lytical expansion at that special point ^U]. A simpler and faster possibility 
(proposed also in is to fit an approximant (the expansion around that 
special point up to the requested precision) for values of p^ where the method 
works, then using it closer to the special point. A test performed in the case 
of the sunrise MI where we can check the results show that the method 
works rather well with some cautions jl3j . 
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Subtracted differential equations are used wlien starting from = co 
or from threshold, as that points are not regular points of the MDE. 

Remarkable self-consistency checks are easily provided by comparing the 
results obtained either starting from the same point and choosing different 
paths to arrive to the same final point, or choosing directly different starting 
points and again arriving to the same final point. 

The execution of the program is rather fast and precise: with an Intel 
Pentium III of 1 GHz we get values with 7 digits requiring times ranging 
from a fraction of a second to 10 seconds of CPU, and with 11 digits from 
few tens of seconds to one hour. 

If A = L/N is the length of one step, L is the length of the whole 
path and the total number of steps, the 4th-order Runge-Kutta formula 
discards terms of order A^, so the whole error behaves as e^K = -/VA^ = 
/N^, and a proper choice of L and N allows the control of the precision. 

Indeed we estimate the relative error, as usual, by comparing a value 
obtained with steps with the one obtained with A^/10 steps, enK = 
\y{N) — V {N /!{})] /V {N) ^ to which we add a cumulative rounding error 
e^^e = y/N X 10"^^ due to our 15 digits double precision FORTRAN im- 
plementation. 

The general massive sunrise MI are numerically well studied in literature 
and several numerical methods are developed, such as multiple expansions 
[Tl| , or numerical integration El ^1 El El Ej ■ Comparisons are pre- 
sented in ^U] with some values available in the literature |141 119j with ex- 
cellent agreement (up to more than 11 digits). For the 4-denominator MI in 
jl3j we obtain complete agreement with and with |2L)lll5j . Calculations 
via numerical integration are also in |21j . 



5. Perspectives 

The presented method for numerically advancing the solutions of the 
MDE is rather precise and competitive with other available methods for 
numerical MI calculations. 

It seems to be possible to complete the 2-loop self-mass for arbitrary 
internal masses and we have completed the 4-denominators case [13j. 

We think that the extension to graphs with more loops or legs do not 
present serious problems, even if the growth in the number of MI increases 
the computing time. 

It is worth to mention that the method relies on the same MDE, which 
are used also for analytic calculations, so it provides a 'low-cost' comforting 
cross-check for those results. 
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